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We present a QCD sum rule analysis for the the newly observed resonance X c (3250) by assuming 
it as a Do(2400)iV molecular state. Technically, contributions of operators up to dimension 12 are 
included in the operator product expansion (OPE). The numerical result is 3.18 ± 0.51 GeV for the 
-Do(2400)7V molecular state, which is in agreement with the experimental data of X c (3250). This 
supports the statement that X c (3250) could be a Dq(2400)N molecular state. The mass for the 
bottom counterpart BqN state is predicted to be 6.50 ± 0.49 GeV. 

I 1 PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 

o: 

> ! I. INTRODUCTION 

O- 

Very recently, BaBar Collaboration reported the measurement of the baryonic B decay B — > 
Y,+ + pTT~Tr~ and observed a new structure in the £+ + 7t~7t~ invariant mass spectrum at 3.25 GeV [l|]. 
\ For simplicity, one could name the new structure as X c (3250). Soon after the experimental observation, 
He et al. have suggested that A" c (3250) could be a Z?q(2400)7V molecular state from an effective Lagrangian 
calculation Q. Theoretically, the molecular concept is well and truly not a new topic but with a history. 
It was put forward nearly 40 years ago in Ref. [3] and was predicted that molecular states have a rich 
spectroscopy in Ref. |4|. The possible deuteron-like two- meson bound states were studied in Ref. @. In 
recent years, some of "X" , "Y" , and "Z" new hadrons are ranked as possible molecular candidates. Such 
as, X(3872) could be a DD* molecular state 0-0; ^(4350) is interpreted as a D* s D* s0 state [ll|, 0; 



YY4260) is proposed to be a XcP° UM 01 an w Xci state [14j; Z + (4430) is deciphered as a D*D\ state 
HEH; #6(10610) and Z b (10650) could be B*B and B*B*, respectively [izlH]- If molecular states can 



> 

r*"** ■ be completely confirmed by experiment, QCD will be further testified and then one will understand the 
QCD low-energy behaviors more deeply. Therefore, it is interesting to study whether the newly observed 
X c (3250) state could be a Dq(2400)N molecular state. Especially, the new resonance may open a new 
window to study molecular states with the meson-baryon configuration. 

In the real world, quarks are confined inside hadrons and the strong interaction dynamics of hadronic 
systems is governed by nonperturbative QCD effect completely. Many questions concerning dynamics of 
the quarks and gluons at large distances remain unanswered or understood only at a qualitative level. It 
is quite difficult to extract hadronic information quantitatively from the basic theory of QCD. The QCD 



sum rule method 19( is a nonperturbative formulation firmly based on the first principle of QCD, which 
has been successfully applied to conventional hadronic systems, i.e. mesons or baryons (for reviews see 
2^] and references therein). For multiquark states, there have appeared fruitful results from QCD sum 



rules these years (for a review on multiquark QCD sum rules one can see [24[ and references therein). In 
particular for hadrons containing five quarks, some authors began to study light pentaquark states in Refs. 
[25j. The application of QCD sum rules to heavy pentaquark states was performed in Ref. [26j for the 
first time. 

In this work, we devote to investigating that whether the newly observed resonance X c (3250) could 
be a Dq(2400)N molecular state in the framework of QCD sum rules. As a byproduct, the mass for the 
bottom counterpart BqN is also predicted. The rest of the paper is organized as three sections. We discuss 
QCD sum rules for molecular states in Sec. [XT] utilizing similar techniques as our previous works [27[. The 
numerical analysis and discussions are presented in Sec. IIII1 and masses oi Dq(2A00)N and BqN molecular 
states are extracted out. The Sec. |IV] includes a brief summary and outlook. 
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II. QCD SUM RULES FOR MOLECULAR STATES 

The QCD sum rules for molecular states are constructed from the two-point correlation function 

n(g 2 ) = z / d 4 xe^(0|T[j(^)J(0)]|0). (1) 



In full theory, the interpolating current for jDq(2400) or Bq meson can be found in Ref. [28[, and the one 
for nucleon N has been listed in Ref. [2^]. One can construct the Dq(2400)-/V or BqN molecular state 
current from meson-baryon type of fields 

3 = (Q c 'q c >)(eabcqT a C^q b 2l , l5 ql), (2) 

where Q is heavy quark c or b, and q\, q2, as well as q% denote light quarks u and/or d. The index T 
means matrix transposition, C is the charge conjugation matrix, with a, 6, c and d are color indices. One 
should note that meson-baryon molecules in the real world are long objects in which the quark pairs are 
far away from each other. The currents in this work and in most of the QCD sum rule works are local and 
the five field operators here act at the same space-time point. It is a limitation inherent in the QCD sum 
rule disposal of the hadrons since the bound states are not point particles in a rigorous manner. 
Lorentz covariance implies that the two-point correlation function in Eq. ([1]) has the form 

ii(q 2 ) = u 1 (q 2 ) + 4iL 2 (q 2 ). (3) 

According to the philosophy of QCD sum rules, the correlator is evaluated in two ways. Phenomenologi- 
cally, the correlator can be expressed as a dispersion integral over a physical spectral function 

^,2-, x2 i + M„ i r°° TmnP hen ( s ) + ^imnP hen ( s ) 

Il(q 2 ) = X H Y72 ^2 + ~ / ds — \ + subtractions, (4) 



m z h - 1 * J so s - i 

where M# is the mass of the hadronic resonance, and Ah gives the coupling of the current to the hadron 
(0|j|ff) = Ajj. In the OPE side, short-distance effects are taken care of by Wilson coefficients, while 
long-distance confinement effects are included as power corrections and parameterized in terms of vacuum 
expectation values of local operators, the so-called condensates. One can write the correlation function in 
the OPE side in terms of a dispersion relation 

n(g 2 ) = r ds^l + nr ic V) + i\ r ds —2 + n 2 0nt V)), (5) 

where the spectral density is given by the imaginary part of the correlation function 

Pi{8) = -ImnP PE ( S ), i = l,2. (6) 

Technically, one works at leading order in a s and considers condensates up to dimension 12. To keep the 
heavy-quark mass finite, one can use the momentum-space expression for the heavy-quark propagator 



^o(p) = "7-^ Tg* A <&(0) , 2 2 , 2 K\{i> + m q ) + {i> + m Q K A 

p — rriQ 4 (p z — tuqY 



- -gH A t B G^{Q)G%(Q) A + _ "5\ B [7atf + mohpW + m Q ) 7/ ^ + m Q ) lv (7) 

4 (p m Ql 

+ ^f ABC G^GlG c e 2 1 + ™ Q )[^(p 2 - 3m 2 ,) + 2m Q (2p 2 - m g )]0 + m Q ). 

48 (p — TUq) 
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The light-quark part of the correlation function can be calculated in the coordinate space, with the light- 
quark propagator 



S a b(x) 



iSab , _ m q 6 ab 



X 2 Sgb 



(gqa ■ Gq) 



T^gG^a^+a^f) 



ix 2 S a b , _ „ . , X 4 8 a b 



12 



48 



-^—^m q {gqa ■ Gq)fi - ^ 10 ^ 3 (qq)(g 2 G 2 ), 



(8) 



which is then Fourier-transformed to the momentum space in D dimension. The resulting light-quark part 
is combined with the heavy-quark part before it is dimensionally regularized at D = 4. Equating the two 
sides for II(g 2 ) and assuming quark-hadron duality yield the sum rules, from which masses of hadrons can 
be determined. After making a Borel transform and transferring the continuum contribution to the OPE 
side, the sum rules can be written as 



Ml /M 2 _ 



d spi ( s )e- s / M2 + sn? ond , 



(9) 



\ 2 H e- M "/ M = 



f° d S p 2 {s)e~ s ' M2 +BTl c 2 ond , 

Jml, 



(10) 



where M 2 indicates the Borel parameter. To eliminate the hadron coupling constant Ah and extract the 
resonance mass Mff, one can take the derivative of Eq. © with respect to 1/M 2 , divide the result by 
itself and deal with Eq. (fTU)) in the same way to get 



M 2 H 



dsp\(s)se 



-s/M 2 



dM-^)Bn? ond ( s )|/|£° d spi {s)e-^ 2 +sn? ond (.s)},(ii) 



M H = 



d SP2 ( S )se-°/ M2 +rf/d(-^)Bnf nd ( S )}/| jf d SP2 {s)e- s ' M2 +5nf nd ( S )},(12) 



where 



/ x pert, x . (qq), \ , (qq) 2 , \ , (gqa-Gq), x . (g 2 G 2 ), x . (g a G 3 ), x . (qq) 3 , x . (qq) (gqa-Gq) , x 

Pi ( S ) = PI ( s ) + Pi ( s ) + Pi ( s ) + Pi ( s ) + Pi ( s ) + Pi ( s ) + Pi ( s ) + Pi ( s ) 



(gqa-Gq) (gqa-Gq) 



(s) + p 



(<i<i)(g 2 G 2 ) 



(s) + p 



(qq){g 3 G 3 ) 



(*) + P 



(g 2 G 2 ) (gqa-Gq) 



(s), 1 = 1,2. 



(13) 



As a matter of fact, many terms of pi(s) are approximate to zero because they are proportional to light 
quarks' masses in the calculations. Thereby, we merely present the spectral densities resulted from Tl2(q 2 ) 
here. Concretely, they can be written as 



pert/ 
P 2 I s 



pT 



P-2 



gqa-Gq) 



pf G \s 
pf G3) (s 



1 



da- g-^— (as — TOq) 4 (q;s + 4mg), 

/A a 

3~2^ A da—^—(as-m Q )(as + m Q ), 

m Q {gqa ■ Gq) f 1 (1 - a) 3 2 2 
2^ A da—^—(as-m Q ) , 

^WG 2 ) f {X-af t ^ 2 
5 • 3 2 • 216^ J A a 



Q"' 



(as - m 2 Q )(as - 2m 2 Q ), 



f^Jl-a) 



-[(as) 2 — 9asmn + 10m 4 
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All)(g<ic-Gq) 

Pi 



(gqa-Gq)(gqa-Gq) 
l J 2 



P-2 



(s 
(s 

P2"" '( s 

(g 2 G 2 )(gqa-Gq) 



(s 



and 



BIl 



cond 



m Q (qqy 



m Q {qq) 
3 • 2 4 7T 2 



3 r 1 



da(l — a), 



(qq) (gqo ■ gg) f 1 , n A2 



2 10 7T 4 



da(l — a), 



m Q (gg)( 5 2 G 2 ) f 1 J (1 - a) 



da- 



3 2 • 2 13 7T 6 

m Q (g<z)(. 9 3 G 3 ) f^ Jl-a) 4 



-[6a 2 (as - mg) + (1 - a) 2 (3as - 4m|)], 



3 • 2 14 7T 6 

TO Q ( g 2 G 2 )(gga-G g ) f 1 (1 - a) 3 
3^2^ ' 



2 6 7T 2 



2 1( V 4 



da— e" m V(^ 2 



n 



m Q (^)( 5 2 G 2 ) 2 ^(1-0)73 



33 . 215^6 



a a 2 M 2 



3 2 • 2 14 7T 6 



e -m^/(aM 2 ) 
2 I r.r.\2 I „2n2\ 



3 3 • 2^ ; "" 



3 3 • 2 10 7T 4 
Q_L-™g/(aM 2 ) 



da- 



* 2 M 2 



m 2 Q /(aM 2 ) 



e 



( g - g )( 5 2 G 2 )(Ma-Gg) r 1 



3 2 • 2!% 4 



da 



a a 2 M 2 

(2 - 4a + 3a 2 )TO| (1 - afm^ 



1 - 



a 3 M 2 



a 4 (M 2 ^ 2 



-m 2 /(aM 2 ) 



* tfG 2 )^ ■ Gq) fi (1-a) 3 



3 2 • 2 13 7T 6 



da- 



. e -™? g /("M 2 ) 



! no,(gq<T-Gq)(g 3 G 3 ) f 1 (1 - a ) 3 f 3 



3 2 • 2 14 7T 6 



a- 



a a 2 M 2 



5_ \„-™ 2 Q /(aM 2 ) 



(14) 



for £)g(2400)7V or Sq^V state. The lower limit of integration is given by A = tuq/s. 



III. NUMERICAL ANALYSIS AND DISCUSSIONS 

In this section, the sum rule (|12p is numerically analyzed. The input values are taken as m c = 1.23 ± 
0.05 GeV, m b = 4.24±0.06 GeV, (qq) = -(0.23±0.03) 3 GeV 3 , (gqa-Gq) = m% (qq), m 2 , = 0.8±0.1 GeV 2 , 
(g 2 G 2 ) = 0.88 GeV 4 , and (g 3 G 3 ) = 0.045 GeV 6 0. In the standard QCD sum rule approach, one can 
analyse the convergence in the OPE side and the pole contribution dominance in the phenomenological 
side to determine the conventional Borel window: on one hand, the lower constraint for M 2 is obtained by 
the consideration that the perturbative contribution should be larger than condensate contributions; on 
the other hand, the upper bound for M 2 is obtained by the restriction that the pole contribution should 
be larger than the continuum state contributions. Meanwhile, the threshold is not arbitrary but 
characterizes the beginning of continuum states. 

On the choice of Borel windows in this work, we would make some more discussions below. It has been 



shown in some Refs., e.g. [30|,|31(, that some condensate contributions are very large, making the standard 
OPE convergence to happen only at very large values of M 2 . Hence one may not find the conventional Borel 
window. In fact, it has appeared the similar problem here. The comparison between pole and continuum 
contributions from sum rule (jlOl) for Dq(2400)N state for ^/so = 3.8 GeV is shown in the left panel of FIG. 
1, and its OPE convergence by comparing the perturbative with other condensate contributions is shown in 
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the right panel. Even if we choose some uncritical convergence criteria, e.g. the perturbative contribution 
should be at least bigger than each condensate contribution, there is no standard OPE convergence up to 
very large values of M 2 . The consequence is that it is unable to find the conventional Borel window where 
both the OPE converges well (i.e. the perturbative contribution bigger than each condensate contribution) 
and the pole dominates over the continuum. 




1.5 2 2.5 3 3.5 2 3 4 5 6 7 8 



FIG. 1: In the left panel, the solid line shows the relative pole contribution (the pole contribution divided by the 
total, pole plus continuum contribution) and the dashed line shows the relative continuum contribution from sum 
rule (|10p for yfso = 3.8 GeV for DJ(2400)iV state. The OPE convergence is shown by comparing the perturbative 
with other condensate contributions from sum rule (|10[) for y/so = 3.8 GeV for Dq(2400)N state in the right panel. 

In order to obtain some useful hadronic information from QCD sum rules, one may take into account 
several possible ways to solve the problem. Firstly, one could try releasing the criterion of pole dominating 
over continuum and take some high values of Borel parameter M 2 . Thus, OPE series may converge well. 
However, there may occur some other problem. As one knows, the phenomenological side of the sum rule 
can be expressed as e.q. (0} due to the "single-pole+continuum states" hypothesis. Once the value of 
of M 2 is chosen too large, the single-pole dominance condition may be spoiled. Secondly, one could try 
pushing the threshold parameter ^fs^ to a very large value, and the maximum value of M 2 will be enhanced 
with the increasing of ^/sq. Thus, one may find the Borel window satisfying both the perturbative bigger 
than condensate contributions and the pole bigger than continuum contributions. However, the threshold 
parameter y/so is not arbitrary but characterizes the beginning of the continuum states. With too large 
values of y/so, contributions from high resonance states and continuum states may be included in the pole 
contribution. Hence, the QCD sum rule may not work normally. Thirdly, one could try releasing the 
strict convergence criterion of perturbative contribution larger than each condensate contribution in some 
case. Here, we consider the ratio of perturbative contribution to the "total OPE contribution" (the sum of 
perturbative and other condensate contributions calculated) but not the ratio of perturbative contribution 
to each condensate contribution. Not too bad, there are denumerable important condensate contributions 
and some of them could cancel with each other to some extent since they have different signs. Most 
of other high dimension condensate contributions are very small. All these factors bring that the OPE 
convergence is still under control at relatively low values of M 2 . What's more important, we find that the 
ratio of perturbative contribution to the "total OPE contribution" does not change much even including 
some condensate contributions higher than dimension 12. From the Borel curves for the £>q(2400)-/V state 
shown in FIG. 2, one can visually see there indeed exist very stable plateaus and one could expect the OPE 
convergence is under control. Thus, we choose some transition range 2.0 ~ 3.0 GeV 2 as a compromise Borel 
window and take the continuum thresholds as y/so — 3.7 ~ 3.9 GeV, and arrive at 3.18 ± 0.41 GeV for 
Dq(2A00)N state. Considering the uncertainty rooting in the variation of quark masses and condensates, 
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4.5 

FIG. 2: The mass of Do(2400)iV state as a function of M 2 from sum rule (j 12|) is shown. The continuum thresholds 
are taken as y^so = 3.7 ~ 3.9 GeV. 
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FIG. 3: The mass of BqN state as a function of M 2 from sum rule (|12[l is shown. The continuum thresholds are 
taken as = 7.0 ~ 7.2 GeV. 

we gain 3. 18 ±0.41 ±0.10 GeV (the first error reflects the uncertainty due to variation of -^io and M 2 , and 
the second error resulted from the variation of QCD parameters) or 3.18 ±0.51 GeV for Z?q(2400)7V state. 
There comes forth the same problem for BqN as the above case for Dq(2400)7V, and we treat it similarly. 
The mass of BqN state as a function of M 2 from sum rule (fT2"j) is shown in FIG. 3. Graphically, one can see 
there have very stable plateaus for Borel curves. We choose a compromise Borel window 4.5 ~ 6.0 GeV 2 
and take y/so~ = 7.0 ~ 7.2 GeV for B^N state. In the work windows, we obtain 6.50 ± 0.29 GeV for B^N 
state. Varying input values of quark masses and condensates, we attain 6.50 ± 0.29 ± 0.20 GeV (the first 
error reflects the uncertainty due to variation of y^so and M 2 , and the second error resulted from the 
variation of QCD parameters) or 6.50 ± 0.49 GeV for BqN state. 

IV. SUMMARY AND OUTLOOK 

Assuming the newly observed structure A c (3250) by BaBar Collaboration as a Dq(2400)-/V molecular 
state, we calculate its mass value in the framework of QCD sum rules. Technically, contributions of 
operators up to dimension 12 are included in the operator product expansion. The final numerical result 
for the Dq(2400)7V state is 3.18 ± 0.51 GeV, which coincides with the experimental value 3.25 GeV. This 
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consolidates the statement that A c (3250) could be a Dq(2400)N molecular state. Additionally, we have 
also studied the bottom counterpart BqN state and predicted its mass to be 6.50 ±0.49 GeV. By analogy 
with Dq(24O0)N state, this bottom counterpart state could be searched in the £{,7r~7r~ invariant mass 
spectrum in future experiments. 
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